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ABSTRACT 

The growing interest in solar twins is motivated by the possibility of comparing them directly to the Sun. To carry on this kind of 
analysis, we need to know their physical characteristics with precision. Our first objective is to use asteroseismology and interfer- 
ometry on the brightest of them: 18 Sco. We observed the star during 12 nights with HARPS for seismology and used the PAVO 
beam-combiner at CHARA for interferometry. An average large frequency separation 134.4 + 0.3 //Hz and angular and linear radiuses 
of 0.6759 + 0.0062 mas and 1.010 + 0.009 R© were estimated. We used these values to derive the mass of the star, 1.02 + 0.03 M©. 

Key words. Stars: individual: 18 Sco - Stars: oscillations - Techniques: radial velocities - Techniques: interferometric - Methods: data 
analysis 



1. Introduction 

Solar twins, defined as spectroscopically identical to the Sun 
( Cay rel de Strobel et al.|198l) , are important beca use they allow 
precise differential analysis relative to the Sun (Ram irez et al. 
2009} |Melendez et al.|2009] ). The brightest solar twin is 18 Sco 
(HD 146233, HIP 79672; V = 5.5), whose mean atmospheric 
parameters are r eff = 5813 + 21 K, logg = 4.45 + 0.02 and 



[Fe/H] = 0.04 + 01 (|Takeda & T ajitsu 2009; Rami rez et al. 
2009 ] |Sousa et al.||2008[ |M elendez & Ramirez 2007 MTakeda 
2006HValenti & Fischer|2005| l. Its ro 



et al. 2007 ; Melendez et al. 



tation rate and magnetic field are also similar to solar ones (Petit 
et al.||2008] ). Its position in the H-R diagram indicates that the 



star should be slightly youn ger and more massive than the Sun 
( |do N asciment o et al.|2009[ and references therein). 

During the past decade, asteroseismology and interferom- 
etry have arisen as powerful techniques for constraining stel - 
lar parameters (e.g., |Cunha et~aL 2007 1 |Creevey et al. 2007). 
Asteroseismology involves measuring the global oscillation 



* Based on observations collected at the European Southern 
Observatory (ID 183.D-0729(A)) and at the CHARA Array, operated 
by Georgia State University. 



modes of a star (which for Sun-like stars are pressure-driven p 
modes). It is the only observational technique that is directly sen- 
sitive to the deeper layers of the stellar interior, since the charac- 
teristics of the modes depend on the regions through which the 
waves travel. Interferometry requires long-baseline interferome- 
ters capable of resolving distant stars, hence allowing measure- 
ment of their radii. 

These techniques have already been combined to study the 
bright sub-giant /3 Hyi ( [North et al.|2007 ), for which a mass was 
derived through homology relations. Here, we apply a similar 
method to 18 Sco. In Section [2] we present the asteroseismic 
data and describe the method used to derive the average large 
frequency separation. In Section [3] we describe the interferomet- 
ric measurements that, combined with the parallax, allow us to 
estimate the radius. In Section[4]we use these quantities to derive 
the mass. 



2. Asteroseismology 

Detecting solar-like oscillations in a fifth-magnitude star from 
the ground is challenging and only a few instruments offer the 
required efficiency and high precision. We observed 18 Sco us- 
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Fig. 1. Time series of radial velocities (upper panel) and their 
uncertainties (lower panel) from HARPS observations of 18 Sco. 



ing the HARPS spectrograph on the 3.6-m telescope at La Silla 
Observatory, Chile (May or et al.|20 03). The data were collected 
over 12 nights from 10 to 21 May 200^] We used the high- 
efficiency mode with an average exposure time of 99.6 s. This 
resulted in a typical signal-to-noise ratio at 550 nm of 158, 
with some exposures reaching as high as 240. The measured 
radial- velocity time series (filtered for low-frequency variations) 
is shown in Fig. [I] We obtained 2833 points with uncertainties 
in general below 2 m s" 1 . The average dispersion per night is 
-1.1 1 m s -1 and can be attributed mostly to p modes. 

The power spectrum, calculated using the measurement un- 
certainties as weights, is shown in Fig. [2j The spectral window 
W, which is the Fourier transform of the observing window, 
w(0 = Z ~ U) (with ti the mid-exposure time of the z-th expo- 
sure), is shown in the inset of Fig. [2] Strong aliases caused by the 
daily gaps appear on both sides of the central peak at multiples of 
±11.57 /iHz. The power spectrum shows a clear excess around 
3 mHz that is characteristic of solar-like oscillations, reaching 
-0.04 m 2 s -2 (corresponding to amplitudes -20 cm s" 1 ). 

The median sampling time was 135.0 s, including the read- 
out (~ 22.6 s). This leads to an equivalent Nyquist frequency of 
3.7 mHz. In Fig.[2j we clearly see a steep rise in power at 7 mHz, 
corresponding to the folded low-frequency increase. This also 
causes the bump that appears between 3.7 mHz and -5.5 mHz, 
which is an alias of the oscillation spectrum. 

The asymptotic relation for high-order p modes is v n j = 
(n + 1/2 + e s )Av, where v n j is the frequency of the mode with 
radial order n and angular degree /, the average large separation 
Ay and a surface phase offset e s (Tassoul 1980). There is a peri- 
odicity of Av/2 in the frequency distribution, implying that there 
will be a local maximum at this value in the autocorrelation func- 
tion (ACF) of the signal y(t), defined by R yy (r) = E\y*(t)y(t + r)] 
(where E is the expectation value and y* the complex conjugate 
of y). The ACF for 18 Sco is shown in Fig. [3] Since the signal is 
irregularly sampled with daily gaps, it was computed by apply- 
ing the Wiener- Khinchine theorem. 

The large separation estimator is thus simply 



Ay: 



2 x [a.rgma.x(R yy * H)] 



(1) 




0.000, 
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Fig. 2. Power spectrum of 18 Sco, evaluated using a weighted 
Lomb-Scargle "periodogram". The vertical grey dashed line 
marks the location of the equivalent Nyquist frequency. The in- 
set shows the spectral window W, normalized to its maximum. 



where 2F is the domain in which we search for this maximum, 
which we set at - [13000,25000] s (i.e. in the frequency 
range 40-80 yuHz), and H the Fourier transform of a filter h(y) 
that truncates the power spectrum. Indeed, only the range 1500- 
3700 yuHz is considered when computing (using an FFT algo- 
rithm) the Fourier transform of the spectrum. We used zero- 
padding to ensure that the ACF was evaluated at points separated 
by an "equivalen t frequency resolution" -0.01 yuHz. As noted by 
Roxbu rgh ( 2009|, the width of h(v) affects the localization of Ay, 
which is a limitation of the method. 

The next step is to obtain information about the statistical 
properties of the estimator Ay that accounts for the noise in the 
data. We can write 



y = x + s, 



(2) 



l An attempt to reduce the daily aliases, with simultaneous ob- 
servations on SOPHIE (Observatoire de Haute-Provence, run ID: 
09A.PNPS.THEA), was unfortunately plagued by bad weather. 



where y = [yo, . . . , y^] are the measured values of the radial ve- 
locity at times to, . . . , x = [xo...,xn] are the true values of 
the radial velocity and e - [sq,..., s n ] is a vector gathering the 
noise contributions from observational errors. Our goal is to es- 
timate the probability density of Ay conditional on y, p(Av\y). 
To do so, we used a Monte Carlo approach to error propaga- 
tion. We assumed that the noise in the data is a series of realiza- 
tions of independent random variables distributed according to 
the Gaussian distributions N = N(x(ti), cr?), with the cr, given 
by the uncertainties on the data. We simulated time series by 
generating new realizations of the noise distributed according to 
the Ni at each t t and adding them to y t . For each artificial set of 
data 3^, we estimated Av". 

Our process for generating the artificial data means that it 
satisfies 

f=y + e'=x + £ + s\ (3) 

with e' the artificially generated noise, and e! and e both real- 
izations of the same distribution at time t t . Ideally, one wishes 
to estimate the large separation from x. One possibility would 
be to estimate it from several measurements, on different tele- 
scopes at the same times t t . A second way would be to generate 
the y from a model reproducing the data y, then perturbing the 
output of this model, rather than the real observations, which is 
the classical procedure of Monte Carlo estimation of parameters. 
Unfortunately, the knowledge of Ay alone does not permit such 
a model to be set up. 
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Fig. 3. Filtered autocorrelation function for the observed data. 
The shaded area marks the interval in which we searched 
for the local maximum corresponding to Av/2. & is the Fourier 
transform. 

It thus has to be assumed that this bias will not be too severe, 
which can be crudely checked graphically with an echelle dia- 
gram (see Fig. [4]). Further estimations of the large separations, 
using individual frequencies, may give us some information on 
its magnitude. However, the error bars are representative of the 
error propagation: were we able to correct for the bias induced 
by random observational noise, we would expect our confidence 
interval on the large separation to be the same. 

Figure[4]shows the results for our Monte Carlo suite of 10000 
time series. It closely follows a Gaussian distribution with pa- 
rameters [A = 134 A //Hz and cr = 0.3 //Hz. The right panel rep- 
resents the echelle diagram of the observations using this value 
for the mean large separation. 



3. Interferometry 

To measure the angular diameter of 18 Sco, which is expected 
to only be about 0.7 mas, we used long-baseline interferome- 
try at visible wavelengths. We used the PAVO beam combiner 
(Prec ision Astronomical Visible Observations; Irela nd et al.| 
12008] ) at the CHARA array (Center for H igh Angular Resolution 
Astronomy; | ten B rummela ar et al.|2 005). We obtained four cal- 
ibrated sets of observations on 18 July 2009 using the S1-W2 
(211 m) baseline. 

PAVO records fringes in 38 wavelength channels centred on 
the R band (A c ^ 700 nm). The raw data were reduced using the 
PAVO data analysis pipeline (M. J. Ireland et al. in preparation). 
To enhance the signal-to-noise ratio, the analysis pipeline can 
average over several wavelength channels, and for 18 Sco we 
found an optimal smoothing width of five channels. Excluding 
four channels on each end because of edge effects, this resulted 
in six independent data points per scan and hence a total of 24 
independent visibility measurements for 18 Sco. 

Table Q] lists the three stars we used to calibrate the visibili- 
ties. We es timated their angular diameters, 0, using the V-K cal- 
ibration of Kervella et al. (2004). Although the internal precision 
of this calibration, as well as the uncertainties in the photometry 
for all three stars, is better than 1%, we a ssume here conservative 
uncertainties of 5% for each calibrator ( [van Belle & van B elle 
2005). These incorporate the unknown orientation and expected 
oblateness in fast rotators (Roy er et al.|2002 ). These were cho- 
sen to be single stars with predicted diameters at least a factor of 
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Fig. 4. Left panel: echelle diagram corresponding to the mean 
value of Av. Right panel: results from the Monte Carlo experi- 
ment. The histogram shows the distribution of the actual realiza- 
tions, and the continuous line the Gaussian with the correspond- 
ing mean and variance. 



two smaller than 18 Sco and to be nearby on the sky (at a sepa- 
ration d < 10°). Each of the four scans of 18 Sco was calibrated, 
using a weighted mean of the calibrators bracketing the scan. All 
scans contributing to a bracket were made within a time interval 
of 15 minutes. The final calibrated squared- visibility measure- 
ments are shown in Fig. [5] as a function of spatial frequency. 

To determine the angular diameter, corrected for limb d ark- 
ening, we fi tted the following model to the data (Hanbury Brown 
|et aL|1974» : 



V 



2 + 3 



(1 



-^)^ + ^/2) 1/2 



.3/2 



(4) 



with x - 7iB6 LD A~ l . Here, V is the visibility, ^ the linear limb- 
darkening coefficient, J n (x) the n-th order Bessel function, B the 
projected baseline, #ld the limb-darkened angular diameter, and 

X the wavelength at which the observations were done. In our 
analysis, we used \ix - 0.607 ± 0.012, which is interpolated at 
the r e ff, logg, and m etallicity of 18 Sco in the R-filter given in 
the catalog of Claret ( |2000| ). For all wavelength channels, we 
assumed an absolute error of 5nm (-0.5%). 

Interferometric measurements are often dominated by sys- 
tematic errors and therefore require a careful analysis of all er- 
ror sources. To arrive at realistic uncertainties for the angular 
diameter, we performed a series of 10 4 Monte Carlo simulations 
as follow. For each simulation, we drew realizations from the 
observed values (assuming they correspond to the parameters of 
Gaussian distributions) for the calibrator angular diameters, limb 
darkening coefficient, and wavelength channels. With these pa- 
rameters we then calibrated the raw visibility measurements and 
fit the angular diameter # LD to the calibrated data using a least- 
squares minimization algorithm. Finally, we generated for each 
simulation a random sample of 200 normally distributed points 
with a mean corresponding to the fitted diameter and a standard 
deviation corresponding to the formal uncertainty (scaled so that 

X I - 1) of the fit. For each Monte-Carlo simulation these 200 
points were stored to make up the final distribution containing 
2 x 10 6 points. This procedure was carried out for all indepen- 
dent measurements in our data. 

The resulting distribution for the diameter of 18 Sco is 
shown in Fig. 5, along with the best-fitting model. The mean 
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Fig. 5. Calibrated squared visibilities for 18 Sco. The red line 
represents the best model. Each symbol type corresponds to 
one scan, the visibility averaged over five wavelength channels. 
Upper-right inset: Same model represented on a larger scale. 
Lower-left inset: distribution for angular diameter (mean and 
standard deviation represented by vertical lines). 

Table 1. Properties of the calibrators used for 18 Sco. 



Star 



3E 



Sp. Type 6 (mas) 



d (deg) 



HD145607 5.415 
HD145788 6.255 
HD147550 6.245 



A4V 
A1V 
B9V 



0.343+0.017 
0.256+0.013 
0.233+0.012 



0.9 
4.2 
6.5 



http : //simbad . u- strasbg . fr/simbad/ 



http : //www . ipac . caltech . edu/2mass/ 



and standard deviation of this distribution yield #ld = 0.6759 ± 
0.0062 m as. Combined with the Hipparcos parallax of 71.94 ± 
0.37 mas ( |van L eeuwen 2007 ), we find the radius of 18 Sco to 
be R/R Q = 1.010 ± 0.009. We conclude that the radius of 18 Sco 
is the same as the Sun, within an uncertainty of 0.9%. 

4. Mass of 18 Sco 



Gough| ( [19 90) pointed out that the homology relation: 

Av oc M 1/2 7T 3/2 , 



(5) 



holds for main-sequence stars even outside the zero-ag e main 
sequence. Other studies have confirmed this picture (e.g., |Stello| 
|et al.|2009| ). 

Applying the ACF method de scribed above to a one -year 
BiSON time series of the Sun ( [Broomhall et al.||2Q09| and 
references therein), we found the solar large separation to be 
135.229 ± 0.003 /jHz. With our radius measurement this gives 
a mass M = 1.02 ± 0.03 M for 18 Sco. The agreement is 
good with the published estimates derived from indirect meth- 
ods, such as comparison between spect ro- or photometric obser- 
vations and stellar evolutionary tracks (|Valenti & Fischer|2005| 
Melendez & Ra mgez||2007l |Takeda et al.||2007[ |Sousa et al.| 
2008Hdo Nascimento et al.|2009| ). 

The assumption of homology is in general well- supported by 
the comparison to models. Considering a small departure from it 
in the form Ay oc (1 + c)M 1/2 ff~ 3/2 (with c« 1, being a function 
of the structure of the star), we then have for a scaling relative 



to the Sun, c = (Av/Av Q )(M/M Q ) l/2 (R/R Q )- 3/2 - 1. For stellar 
ages characteristics of those quoted for 18 Sco (the dependence 
on the mass and the metallicity being weak), this quantity may 
contribute to an additional ~0.2%-0.4% on the total error on the 
mass. 

The impact of filtering the ACF is not completely negligi- 
ble, and if we vary the lower limit of between 1500 yuHz and 
2000 //Hz, the final estimate may vary by -0.01 M , emphasiz- 
ing the need for individual frequency measurements. 

5. Conclusion 



We presented the first asteroseismic and interferometric mea- 
surements for the solar twin 18 Sco. These allowed us to esti- 
mate a mass for this star independent of the previous spectro- 
photometric studies, which are still being confirmed. This work 
shows the possibilities offered by asteroseismology, even from a 
ground-based single site, and iby nterferometry. Our results con- 
firm that 18 Sco is remarkably similar to the Sun in both radius 
and mass. 

The next step will involve measuring the individual oscil- 
lation frequencies and performing full modelling using all the 
available observations. It will hopefully reduce the uncertainty 
on the estimated age, improving our k nowledge of the physical 
state of 18 Sco (do Nascime nto et al.||2009| ). This will provide 
a more precise picture of its interior and give information on 
the depth of its external convective zone ( Monteiro et al.|2000 ), 
which is necessary if o ne wishes to study its magnetic activity 
cycle ( [Petit et al.|2008] ). 
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